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We study the phase diagram of strongly interacting matter with three quark flavors at 
■ low and intermediate densities and non-zero temperatures in the framework of an NJL-type 

' model with four-point interactions. At large densities, when the interactions are weak due 

to asymptotic freedom, quarks form loosely bound Cooper pairs. However, when the density 
decreases, interactions become stronger and quark Cooper pairs transform smoothly into 
tightly bound diquark molecules. We flnd that such molecules are stable at low density 
and temperature and that they dissociate above a temperature T^iss of the order of the 
chiral phase transition temperature Tc ~ 170 MeV. We also explore the conditions under 
which these molecules undergo Bose-Einstein condensation (BEC). We find that BEC is 
only possible if we increase the attractive interaction in the diquark channel to (probably 
unrealistically) large values. 



^ ' §1. Introduction 

\D 

At asymptotically large densities, due to asymptotic freedom of QCD/) the in- 

■ teraction between quarks becomes arbitrarily weak. Then, single-gluon exchange 
. constitutes the dominant contribution to the quark-quark interaction. Single-gluon 
I exchange is attractive in the color-antitriplet channel. Thus, we expect the forma- 

■ tion of quark Cooper pairs which destabilize the Fermi surface^-* and lead to color 
superconductivity,'^) see^^ for reviews. 

At sufficiently large densities, the coupling constant is sufficiently small so that 
^ ' color superconductivity can be analyzed rigorously in the framework of QCD, using 

resummation techniques based on perturbative methods and power counting. 
This treatment is the analogue of weak-coupling BCS theory in condensed matter 
physics. The zero-temperature gap parameter turns out to be parametrically small 
in the coupling constant, (/>o ~ /i exp(— where ^ is the quark chemical potential. 
This, in turn, leads to "large" Cooper pairs, i.e., the quark correlation length ^ ~ 
^ ~ fi^^exp{l/g) is parametrically larger than the interparticle distance ~ /x~^.^) 
When the density decreases, however, the strength of the quark-quark interaction 
and, thus, the gap parameter, increase, and the correlation length decreases. Let us 
note that a signature of strong correlations in the normal phase is the appearance 
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of a so-called pseudogap in the vicinity of the transition to the superconducting 
phase.^^^'^-*^) When the correlation length ^ < quark Cooper pairs should be 
regarded as tightly bound diquark molecules. These molecules are stable at all 
temperatures below their so-called dissociation temperature T^iss-^^-* Above T^iss, 
they decay into the two quarks constituting the molecule. At small temperatures, 
diquark molecules may undergo Bose-Einstein condensation (BEC).^)'^^)"^'^^ 

In these proceedings, we investigate bound diquark states and the possibility 
that they undergo BEC in the phase diagram of the quark matter. We use an NJL- 
type model with four-quark interactions. The strength of the attractive diquark 
interaction in the color-antitriplet channel is regarded as a free parameter. We show 
that bound diquark molecules appear at low densities and temperatures for all values 
of the diquark coupling strength studied here. We also find that BEC of diquarks 
occurs for (probably unrealistically) large values of the diquark coupling. 

The remainder of this work is organized as follows. In Sec. 2 we introduce the 
model and the formalism in order to study diquark correlations. In Sec. 3 we present 
our numerical results. We conclude our work in Sec. 4. Our units are h = = c = 1, 
the metric tensor is g'^'^ = diag(+, —,—,—). 

§2. Model and Formalism 

In this work, we employ an NJL-type model Lagrangian for three quark flavors, 



a=0 



7,c 



where the quark field tp'^ has color, a = r,g, b, and flavor, a = u,d, s, indices. The 
current quark mass matrix is given by m = diagj(m„, m^, ms), and are (twice) 
the generators ofU{3). The charge-conjugate spinors are V'c = C'ljj'^ and = i^'^C, 
where C = 27^7° is the charge conjugation matrix. A sum over doubly appearing 
upper and lower indices in color and flavor space is implied (but not if both are 
upper or lower indices). 

In the Lagrangian (2-1), the terms proportional to Gs are quark-antiquark four- 
point interactions in the scalar and pseudoscalar channel, respectively. The terms 
proportional to Go parametrize the diquark four-point interaction in the color- 
antitriplet, flavor-antitriplet channel. For one-gluon exchange in QCD, this channel 
is attractive and thus leads to color superconductivity. Note that the diquark term 
can also be obtained from the quark-antiquark term by a Fierz transformation; in 
this case, the diquark coupling strength is fixed, Gd = 0.75Gs-'^^^ For the sake of 
simplicity, we neglect the eflPect of the U(1)a anomaly, so there is no t'Hooft-type 
six-point interaction term in Eq. (2-1). 
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In mean-field approximation, the thermodynamic potential is 



EiG^ + i: 8Gs -2^:7 (^^A/,cln[^ (^a.„,p)], (2-2) 



c=l a=l 

where 



Ma = ma-^Gs{lljatpa), (2-3) 

/Ae = 2Gz.(C(Pctf V'J), (2-4) 

are the constituent quark masses and the gap parameters for color superconductivity, 
respectively, with (-Pc)af = ijs^'^^ c^abc- The 72 x 72 Nambu-Gor'kov propagator is 
defined by 

with ^ = iujn'jo ~ P ■ 7- Here, a;„ = (2n + l)TrT, n = 0, ±1, ±2, . . ., arc the fcrmionic 
Matsubara frequencies. In this first exploratory study of diquark molecules and BEC, 
we assume the quark chemical potential /Lt to have a common value for all flavors. In 
application to compact stellar objects, this changes due to the conditions of overall 
electric (and color) neutrality, as well as f3 equilibrium. In application to heavy-ion 
collisions, this is also diflFerent, as (for isospin-symmetric nuclei) fJ-u = fJ'd = fJ-i while 
lis = 0. 

The physical values of the variational parameters Ac and Mq, satisfy the gap 
equations 

— - = and — — = 0. 2-6 

dAc dMa ^ ' 

In mean-field approximation, the phase transition from the color-superconducting 
to the normal phase is of second order. This means that the order parameter for 
condensation, A^ goes to zero smoothly as a function of temperature and the critical 
temperature for the transition can be determined from 

' 0. (2-7) 



Ac OA, 

Since up and down fiavors are treated symmetrically in our model, we assume 
M„ = Md and Ai = A2. Because of explicit chiral symmetry breaking by a nonzero 
current quark mass, {tpa'ipa) is nonzero for all T and fi. In the following, we refer to 
the phase with A^ ^ and Ai^2 = as the 2SC phase, and As ^ and Ai^2 7^ 
as the CFL phase.^^ The phase with Ai = A2 = = corresponds to unpaired 
quark matter. 

At nonzero temperature, the order parameters Ac and fluctuate around their 
mean-fleld values. The propagation of these fluctuations in unpaired quark matter 
is characterized by the retarded propagator 

(x, t; x', t') = -ie{t - t')([V^(x, t)PcV'c(x, t), ^c(x', t')Pc^{^', t')]) 

^"^'^ D^iLU, k) e--(*-*')+*k.(x-x')^ (2-8) 



/ 



(27r)^ 
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where c = 1, 2, or 3 denotes the down-strange, up-strange, or up-down diquark field, 
respectively. In the random phase approximation, the diquark propagators are given 

by 

1 Q^(^,p) 



2 l + GDQ§{oo,py 



(2-9) 



where Q^{lo, p) is the one- loop quark-quark polarization function. For imaginary 
energies u = ivn, where f„ = 2mrT, n = 0, ±1, ±2, 
frequencies, it is given by 



are the bosonic Matsubara 



Qc{iiyn,p) = 2T 



y f— 



e/^|TrD,c[^/3(«WTO, ({)Qi{wn + ii^m, P + q)], (2-10) 



with the trace taken over Dirac and color indices. Here, Qa{i^ni p) = [{i^n + A*)7o ~ 
p • 7 — Ma]~^ are the Matsubara Green's functions for quarks of flavor a. Substi- 
tuting these Green's function into Eq. (2-10) and taking the analytic continuation 
Qcil^^P) = Qc{i^n,p)\iv^-^Lo+ir„ wc obtain 



/3,7 



(2^ 



s,t=± 



{sEp + tE^f - IpP - 5Ml 
EpE^ 

f{tE^-,i)-f{-sEf^ + n) 
Lo + 2fx — sEp — tE^ + iri ' 



(2-11) 



where E 



/3 



q-p|2 + M|, E 



q|2 + M2, 5Mc = \Mf3 - and f{E) 



[exp(£'/T) -|- is the Fermi-Dirac distribution function. The imaginary part of 
Q^{u>,p) denotes the diff'erence of decay and production rates of the diquark field. 
At p = 0, it is given by 



ImQ,«(u;,0) 



13,1 •> y > 



5Ml 



EpEj 



S{u; + 2/1- Eg - E^ 



+ [(1 - /^-)(1 - /-) - /^-/-] Sico + 2i^ + Ef, + E^) 
/^(l - /+) - (1 - /^-)/+] 5{u + 2^i + Ep- E^) 
/+(! - f-) - (1 - /+)/-] 5{u + 2pi-Ef, + E^)] , (2-12) 

where = {exp[(E'Q, =F /u)/T] -|- The first (second) term in curly brackets 

corresponds to the decay of a diquark into two quarks (anti-quarks) and assumes 
nonzero values for uj > 2Mc - 2/x and co < -2Mc - 2/x, with Mc = {Mp + M^) /2. The 
third and fourth terms represent Landau damping of a diquark. They are nonzero 
for —6Mc — 2n<u< 6Mc — 2ji. The numerical results show that the latter processes 
do not affect the stability of diquark excitations in our model, since the pole of the 
diquark field never appears around energies where Landau damping occurs. 
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The poles of the diquark propagator are determined by solving D^{uj, p) = 
or, equivalently, 

1 + GdQ^{uj,p)=0. (2-13) 

Setting a; = |p| = in this equation, one can show that Eq. (2-13) is equivalent to the 
condition Eq. (2-7). This means that there exists a pole of the diquark propagator 
for zero energy and momentum at the critical temperature Tc of a second-order phase 
transition. This fact is known as the Thouless criterion.^^^ Above Tc, the pole moves 
continuously to the fourth quadrant of the complex w-plane. The corresponding 
mode is called the soft mode. If Mc < /x at T = Tc, = is in the continuum and 
the soft mode has a decay width above Tc.^^^ On the other hand, if Mc > 1^, the soft 
mode does not have a decay width and the pole stays on the real axis. Then, this 
mode is nothing but a stable diquark molecule. -^^^ As T increases, the pole will move 
along the real axis until, at the dissociation temperature T^jgg, it eventually arrives 
at the threshold of the decay process into two quarks u^^^ = 2(Mc — fj,) . (Note that 
Mc — /X is the energy required to put one additional quark into the system. If the 
diquark energy exceeds twice this value, the diquark will decay.) The dissociation 
temperature T^j^^ is determined by solving 

1 + GoQci'^Mc - 2/x, 0) = 0. (2-14) 

If stable diquarks are formed above Tc, it is natural to associate the superfluid phase 
below Tc with BEC of these diquark molecules. In the following, we regard the 
region of a color-superconducting phase where fi < Mc is satisfied as a Bose-Einstein 
condensed phase. ^^-^ Note, however, that this is just a rough estimate to separate 
BEC and BCS regimes; these two limits are connected continuously and there is no 
sharp phase boundary. 

§3. Numerical Results 

In this section, we present the phase diagram in the T-fi plane for several values 
of the diquark coupling constant Gd- In Fig. 1, we show the phase diagram for 
Gd/Gs = 0.75, which is the canonical value arising from a Fierz transformation of 
the NJL-type four-point interaction.^^)' One observes two color-superconducting 
phases, the 2SC and CFL phases, at high ^ and low T. At T = 0, these phases 
are separated by a first-order phase transition; the first-order transition terminates 
at nonzero temperature. Also shown in Fig. 1 are the dissociation temperatures for 
stable diquark molecules; below the curve labelled T^^^^, bound states of up and down 

1 2 

quarks are stable, and below T^^^^, bound states of up with strange and down with 
strange quarks are stable. 

As discussed in the previous section, BEC of diquark molecules requires that 
H < Mc inside a color-superconducting phase. In Fig. 1, we show where ji = Mc, the 
regions to the left of these lines satisfy ^ < Mc- One sees that these lines terminate 
at the first-order transition between normal and superconducting phase and that 
fi < Mc is not satisfied inside a color-superconducting phase. Therefore, there is no 
BEC for Gd/Gs = 0.75. 
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Fig. 1. The phase diagram in the T-/i plane for Gd/Gs ~ 0.75. Bold and thin solid lines represent 
first- and second-order phase transitions, respectively, while dotted lines denote the dissociation 
temperature of diquark molecules formed of up and down quarks, T^i^^, and of up and strange 
as well as down and strange quarks, T^i^s' respectively. The condition fj, = Mc is denoted by 
dash-dotted lines. 



In Fig. 2, we show the phase diagram for Gd/Gs = l-l- We see that the 
regions of the 2SC and CFL phases expand towards lower /x and higher T. Now 
there appears BEC inside the 2SC phase, namely where < M3, shown by the 
shaded area in Fig. 2. One also sees that the dissociation temperatures T^^^^ increase 
with Gd- The values of T^^^^ at fj, = are comparable to or even larger than the 
critical temperature of the QCD phase transition as determined by lattice QCD 
calculations, Tc^lqcd — 170 MeV.^^^ Therefore, bound diquarks can survive even in 
the quark-gluon plasma phase if the diquark coupling is sufficiently large. 

The other interesting feature shown in Fig. 2 is the fate of the first-order phase 
transition. One sees that the corresponding transition line terminates at smaller tem- 
peratures for Gd/Gs = 1-1 than for Gd/Gs = 0.75. This result can be interpreted 
as a result of the interplay between chiral symmetry breaking and color superconduc- 
tivity. As far as we have checked, the endpoint of the first-order transition smoothly 
approaches T = {) as Gd increases. *^ 

Finally, let us investigate the phase diagram for an unrealistically large diquark 
coupling of Gd/Gs = 1-5, cf. Fig. 3. We observe that the BEC region becomes 
rather wide. The critical temperatures Tc for the 2SC and CFL phases and T^^^^ also 
increase substantially. If the diquark coupling is raised further, even the vacuum at 
T = fi = becomes a Bose-Einstein condensate of diquark molecules. 



We do not observe a new endpoint at a lower temperature, found in Ref. 
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Fig. 2. The phase diagram in the T-^ plane for Gd/Gs = l-l- 
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Fig. 3. The phase diagram in the T-fi plane for Gd/Gs = 1-5. 



§4. Conclusions 

In this paper, we studied the phase diagram of three-flavor quark matter with 
particular emphasis on the formation of bound diquark molecules and the possibil- 
ity that they undergo Bose-Einstein consensation (BEC). We found that diquark 
molecules appear at low densities and temperatures smaller than their dissociation 
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temperature Tdiss- This dissociation temperature may exceed the deconfinement tem- 
perature and imphes that diquarks can be tightly bound even in the quark-gluon- 
plasma phase. We also found that BEC of diquarks is possible for large diquark 
interaction strengths Gd- This finding is in agreement with the results of Ref.-*^^) 
However, the required values for Go are probably unrealistically large, which makes 
it rather unlikely that a BEC of diquarks will be observed in compact stellar objects 
or even heavy- ion collisions. 
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